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1. Introduction 

Statistical classification of remotely sensed images attempts to discriminate between surface 
cover types on the basis of the spectral response recorded by a sensor. It is well known 1-3 that 
surfaces reflect incident radiation as a function of wavelength producing a spectral signature specific 
to the material under investigation. Multispectral and hyperspectral sensors sample the spectral 
response over tens and even hundreds of wavelength bands to capture the variation of spectral re- 
sponse with wavelength. Classification algorithms then exploit these differences in spectral response 
to distinguish between materials of interest. Sensors of this type, however, collect detailed spectral 
information from one direction (usually nadir); consequently, do not consider the directional nature 
of reflectance potentially detectable at different sensor view angles. 

Improvements in sensor technology have resulted in remote sensing platforms capable of detect- 
ing reflected energy across wavelengths (spectral signatures) and from multiple view angles (angular 
signatures) in the fore and aft directions. Sensors of this type include: the moderate resolution 
imaging spectroradiometer 4 (MOD IS ), the multiangle imaging spectroradiometer 5 ( MIS R), and the 
airborne solid-state array spectroradiometer 6 (AS AS). 

A goal of this paper, then, is to explore the utility of BRDF models in the selection of optimal 
view angles for the classification of remotely sensed images by employing a strategy of searching for 
the maximum difference between surface BRDFs. After a brief discussion of directional reflect ante 
in Section 2, attention is directed to the Beard-Maxwell BRDF model and its use in predicting 
the bidirectional reflectance of a surface. The selection of optimal viewing angles is addressed in 
Section 3, followed by conclusions and future work in Section 4. 

1.1. NEFDS Spectral Database 

A collection of spectral datasets was obtained from the National Imagery Resource Library 
(NIRL) and used in this study. Intimately related to the Materials Exploitation Database (MED), 
the Spectrum Archival Library (SAL), and the Spectral Catalog, 7 the Nonconventional Exploitation 
Fact ors Data System (NEFDS ) database is composed of spectral reflect ante measurement data in 
the visibl^near-infrared (VNIR), the mid-infrared (MIR) and combined visible/near-mid-infrared 
(VNMIR) for samples of selected materials. 

Spectral samples of a concrete runway, a galvanized steel rooftop, and a couple of painted 
surfaces in the visible/near-infrared (VNIR) were selected for input to a BRDF model for the 
pairwise estimation of optimal view angles in the classification of remotely sensed images. Detailed 
properties of the materials used in this analysis have been assembled and are presented in the 
Spectral Catalog cited earlier. 

*This work was supported in part by NASA ACE under contract #NCCW-0087. 
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2 , Directional Reflectance 


The directional-hemispherical reflect ante defined as 


Pdh 


d$ r 


( 1 ) 


is the fraction of incident radiant flux density from direction, that is reflected by the surface 

into all possible directions of the hemisphere surrounding the sample. By definition, then, p' dh 
does not provide any information on the directional nature of reflectance and is only adequate for 
surfaces that reflect isotropically. 


Many surface materials of interest, however, are anisotropic reflectors so that the spectral re- 
sponse is also dependent on the viewing geometry of the sensor. 8-10 This fact has lead to classifica- 
tion difficulties in which single view imagery was used to discriminate between surface materials . 11 
A more general and useful description for the directional nature of reflected radiant flux is con- 
tained in the bidirectional reflectance distribution function (BRDF) which relates the directional 
distribution of exitant flux to the incident radiant flux striking the surface. 


2 . 1 . Bidirectional Reflectance 



Figure 1 . Co-ordinate System Defining Solar Angles and Sensor Viewing Angles. 


The BRDF is a function of four angles (Figure 1) and is defined as the ratio of reflected 
radiance from a surface to the irradiance incident to the surface from an illuminating source. This 
relationship is seen in Equation 2 


Pbd(®i> ©r! 


dL r (Qj , Q r ; A) 
dEi{Qi; A) 


( 2 ) 


where the dependence on incident angles, ©;={0;,<fo}, exit ant angles, 0 r = f9„ cf> T }, and wave- 
length, A, is shown. 
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2 . 2 . BRDF Models 


Reflection of light from a surface is a complex phenomenon which can be very difficult to 
explain and accurately predict. Even with a firm grasp of the physics involved, the dynamics of a 
changing environment complicates matters so that a complete understanding is often not possible. 
As a result, several types of BRDF models have emerged; those based on first principles such as 
radiative transfer theory, 12-14 geometrical optics , 15 - 16 physical optics (wave theory of light) and 
those that fit analytic equations to the observed reflectance. 17 

Two mechanisms of reflect ante are generally identified: surface and subsurface (volumetric) 
reflect ante. Surface reflect ante occurs at the interface between materials; usually the air-material 
interface and can be quite anisotropic. Subsurface reflectance is typically isotropic and occurs when 
incident light penetrates the surface and collides with inhomogeneities (such as paint pigments) 
suspended in the substrate. For rough matte surfaces, the observed reflectance is usually diffuse 
(Lambertian) and therefore nondirectional. The BRDF in this case is independent of both the 
incident and exitant directions and can be expressed as 

®r! A) = /> 6d (A). (3) 

Perfectly smooth planar surfaces on the other hand behave as specular reflectors that reflect light 
in the critical mirror angle only. In this instance, the BRDF is highly dependent on the incoming 
and outgoing directions and can be expressed as a Dirac distribution. 18 Most practical surface 
mat erials are somewhere between the two extremes of diffuse and specular reflectors. Therefore, 
BRDF models typically represent the surface reflectance in terms of both a diffuse component and 
a specular component. 

The aim of BRDF models, then, is to predict the reflectance behavior of surfaces under condi- 
tions of incident radiant energy striking the surface from various directions. One such model is the 
Beard-Maxwell BRDF model. 


2,2.1. Beard- Maxwell BRDF 

The Beard-Maxwell (B-M) BRDF model 19 was originally developed to characterize the re- 
flectance properties of painted surfaces but has been used successfully in estimating the reflectance 
of other surfaces as well. Empirically based, this model has seven input parameters that are derived 
from a series of reflectance measurements collected under the controlled conditions of a laboratory 
setting. 

A functional description of the B-M BRDF model is specified by Equation 4 


Pbdi®i> ®r) — 


R(P) Pf a cos 2 0N 1 ~F On 

R( 0) COS0iCOs9 T 1 -f ^SL e -2^/r 


+ Pd + 


2 Pv 

cos8i + cos6 r 


(4) 


where the first term in Equation 4 is the first surface reflectance attenuated by the shadowing and 
obscuration function enclosed in square brackets. The second term, p d , is the diffuse (Lambertian) 
component and the third term is the component due to subsurface volumetric scattering. 


3. Selection of Optimal Viewing Angles 

In this section, the utility of BRDF models in the selection of optimal viewing angles for 
classification of remotely sensed images is explored. To motivate this discussion, consider the B-M 
BRDF of a galvanized steel rootop surface generated under the following conditions: the material 
is oriented in a horizontally flat position ( 0 mati = 00), the solar position is specified by a zenith 


755 




(a) p'bd (b) p's 


Figure 2. B-M BRDF:0526stla (galvanized steel rooftop) 

angle of 6 % =30° and an azimuth of fa =0°, with the BRDF evaluated at a spectral wavelength of 
A=0.635^m. Note that the observed BRDF seen in Figure 2(a) is dependent on viewing angles 
in both the azimuth and zenith directions due to forward scattering (sharp peak on the left) and 
backward scattering (peak on the right) effects. These effects are also seen in the surface component 
of the B-M BRDF shown in Figure 2(b) where once again the forward scattering peak is sharp and 
to the left with a significant backscattering lobe directed toward the solar position. Analagous 
to the variation in spectral signatures among material types, the spatial distribution of reflectance 
with angle gives rise to angular signatures specific to material type. Differences in angular signature 
can lead to an increased ability to distinguish between materials if identified and incorporated into a 
classification system. Our goal, then, is to identify regions of maximum difference between angular 
signatures through the use of the B-M BRDF model. 

The B-M BRDF of several surfaces was generated with the the material positioned horizontally 
flat, the solar position given by a zenith of = 30° and an azimuth of fa = O"; wavelength 
ranged from 0.3- 1.0 /im in increments of 0.005 fim. Positioning of the sensor varied from fa =0" 
to fa = 180° in the azimuthal direction at 5° increments and in the zenith from 0 T =0° to 
& r = 65° every 5° corresponding to each azimuthal direction. As an illustration, B-M BRDF 
surfaces of a galvanized steel rooftop with the sensor positioned at an azimuth of fa = 180° and 

fa — 175° are shown in Figures 3(a) and 3(b), respectively. Under similar conditions, B-M BRDF 
surfaces were generated for an aluminum painted low emissivit y green surface ( 0537 UUUPNT), a 
gray unweathered polyurethane paint on aircraft surface (0741 UUUPNT), and a concrete runway 
surface (0671 UUUCNC). 

Pairwise separability was then computed on the basis of a pointwise root-squared (rs) difference 
between B-M BRDF surfaces of each material. Therefore, corresponding to each azimuth direction 
a root-squared difference surface was determined as shown in Figure 4(a). Inspection of the rs 
difference surface indicates that the maximum difference occurs in the zenith direction of 9 r =30°. 
A cross-section of the surface along 0 T =30° (Figure 4(b)) further shows that the maximum rs 
difference is achieved in the shorter wavelengths around 0.4^m and decreases as wavelength gets 

longer. A similar analysis in each azimuthal direction, fa = {O, 5,10 60,65}, resulted in the 

maximum difference between BRDF surfaces occuringat fa =0°, 9 r = 30°, and in the wavelength 
range of 0.4 - 0.5 jim. 

A separability analysis to determine the optimal viewing conditions for discrimination between 
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(a) <j> T = 180° (b) <j> T — 175° 


Figure 3. B-MBRDF Surface: 0526UUUSTLa (galvanized steel rooftop) 


o.oe. 



(a )( A r = 180° (b) ij) T = 180°, 8 r = 30° 


Figure 4. Pointwise Difference Surface: 0526 UUUSTLa (galvanized steel rooftop) and 0537UU- 
UPNT (green painted aluminum surface) 

a concrete runway (0671 TJUUCNC) and a painted aircraft surface (0741 UUTJPNT) was also con- 
ducted. Results indicated that the maximum rs difference occurred in the primary backscatter 
direction, cp r =0° and optimal 6 r in the range 20° - 40°. Additional analysis indicated that the 
optimal wavelength range for discrimination, however, was in the longer wavelengths from approx- 
imately 0.8 n m to 1.0 pm. 


4. Conclusions and Future Work 

The utility of BRDF models for pairwise discrimination between surface materials through the 
selection of optimal view angles was explored. Employing a strategy of maximum rs difference 
between the respective BRDFs to indicate maximum separability, several cases were considered of 
which two were reported in this paper. I n every case studied, the maximum rs difference occurred in 
the principal plane of the Sun with the BRDF due to backscatter effects yielding a slightly greater 
separability than the forward scattering direction. 


5 



It should be noted, however, that there were several limitations in this preliminary study: our 
analysis at this point was qualitative, atmospheric effects were completely ignored (we expect that 
these effects are not significant at low altitues but could be extremely important at high altitudes), 
and only one solar position was considered. 

Future efforts will concentrate on quantifying the search for optimal viewing angles, incorporat- 
ing atmospheric effects, and exploring the affect that solar position has on the separability y between 
surface BRDFs. 
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